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Abstract. - We examine vortex flow states in periodic square pinning arrays with one row and 
one column of pinning sites removed to create an easy flow crossed channel geometry. When a 
drive is simultaneously applied along both major symmetry axes of the pinning array such that 
vortices move in both channels, a series of coherent flow states develop in the channel intersection 
at rational ratios of the drive components in each symmetry direction when the vortices can 
cross the intersection without local collisions. The coherent flow states are correlated with a 
series of anomalies in the velocity force curves, and in some cases can produce negative differential 
conductivity. The same general behavior could also be realized in other systems including colloids, 
particle traffic in microfluidic devices, or Wigner crystals in crossed one-dimensional channels. 
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I A wide variety of different types of vortex commensu- 
rability and dynamics can be realized in superconductors 
containing a periodic array of artificial pinning sites [1,2]. 
For fields at which the number of vortices is an integer mul- 
tiple or rational fraction of the number of pinning sites, 
various types of vortex crystalline states occur which are 
associated with peaks in the critical current [1-7]. When 
there are more vortices than pinning sites, it is possible to 
have highly mobile interstitial vortices between vortices 
located at the pinning sites [2,5,8,9]. One-dimensional in- 
terstitial vortex motion between pinned vortices has been 
realized in systems with artificially fabricated weak pin- 
liing channels [10-13]. In these channel geometries, oscil- 
lations in the critical current and resistance can arise as a 
function of vortex density due to changes in the number 
of vortex rows moving within the channels. 

One aspect of vortex dynamics in periodic pinning ar- 
rays that has not been studied is vortex motion in systems 
where groups of pins are removed to form easy flow chan- 
nels for interstitial vortices. In this work we examine the 
vortex motion in a system with a square pinning array 
where a row and a column of pins are removed to create 
intersecting channels. Vortices moving in the two channels 
are forced to interact at the channel intersection and can 
experience interference phenomena. Such a pinning ge- 
ometry is very feasible to create experimentally, and the 
removal of individual pinning sites to create diluted pin- 
ning arrays has already been achieved [14]. Driving of 
vortices with crossed external drives and measurement of 
the perpendicular voltage responses has also been demon- 



strated [15]. 

We first study the effective matching field where the 
number of vortices equals the number of pinning sites in 
the unmodified square array and there are an equal num- 
ber of interstitial vortices in each channel. We fix an ex- 
ternal drive along one of the channel directions such that 
the vortices in that channel flow, and then slowly increase 
a second drive in the perpendicular direction, causing the 
vortices in the second channel to flow and interact with 
the moving vortices in the flrst channel. We flnd a series 
of transitions between disordered flows where the cross- 
ing vortices collide and ordered or coherent flows where 
the vortex motion through the intersection is synchronized 
and no collisions occur. The ordered phases create a series 
of anomalies in both velocity components. We also flnd 
negative differential resistance upon passing into and out 
of the ordered phases when vortices in one channel must 
slow down as the drive is increased in order to maintain an 
ordered flow. The overall structure of the transport curve 
is a devil's staircase. It is distinct from the devil's stair- 
case structures found for vortices [16] or colloids [17, 18] 
moving over periodic pinning arrays. In Refs. [16-18], the 
velocity vector locks to symmetry directions of the pinning 
array, so the devil's staircase structure occurs even in the 
limit of a single isolated particle moving over the array. 
In our system, the locking phases vanish in the limit of 
a single interstitial particle moving in the channels and 
are replaced by what we term a vortex gating effect in 
which the motion abruptly switches from one direction to 
the other. The dynamics of vortices in structured pin- 
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Fig. 1: Pinning site locations (open circles), vortex locations 
(filled dots), and vortex trajectories (lines) at B = for a 
square pinning array with one row and one column of pins re- 
moved. The vortex density matches the original square pinning 
array density B^. A constant drive of = 0.2 is applied in 
the y direction and an increasing drive is applied in the 
X direction, (a) At F^P = 0.025, the vortices in the x chan- 
nel are pinned, (b) At F^ = 0.04, the vortex trajectories 
at the intersection are disordered, (c) At F^ = 0.106 an or- 
dered or coherent flow phase occurs, (d) Disordered flow at 
= 0.1105. (e) Ordered flow at = 0.2. (f) Disordered 
flow at Fi* = 0.224. 



ning geometries can also be generalized to other systems 
including colloids in periodic trap arrays [17, 19] or nar- 
row channels [20], particle or bubble flow in microfluidic 
devices [21], and Wigner crystals in narrow channels [22]. 
We expect that the results we describe here can also be 
realized in these other systems where intersecting quasi- 
one-dimensional particle states can be constructed. 

Simulation- We consider a two-dimensional system with 
periodic boundary conditions in the x and y directions. 
The sample geometry is illustrated in fig. 1(a). A single 
column and a single row of pinning sites are removed from 
a square pinning array to create an intersection at the 



center of the sample. The number of vortices N.^, is equal to 
the number of pinning sites that would have been present 
in the original square pinning array, and is larger than the 
actual number of pinning sites Np in the modified square 
array. That is, the applied field B = B^, where B^ is the 
matching field of the undiluted square array. All of the 
pinning sites arc occupied and the remaining vortices are 
located in the interstitial regions in the empty row and 
column. A single vortex i located at position obeys 
the overdampcd equation of motion: 



'dT 



(1) 



The repulsive vortex-vortex interaction force is given by 
F™ = J2f^t^Ki{R,j/X)fi,j, where Ki is the modi- 
fied Bessel function, A = 0Q/(27r/xoA'^), 0o ~ h/2e is 
the flux quantum, A is the London penetration depth, 
Rij = [Ri — Rj l, and Ry = (Ri — Rj)/i?y . The damping 
constant rj = (/)§c?/27r^^/9jv, where ^ is the superconduct- 
ing coherence length, pat is the normal state resistivity 
of the material, and d is the thickness of the supercon- 
ducting crystal. The pinning force arises from attrac- 
tive parabolic potential wells arranged as shown in fig. 1 
with a maximum force of Fp and a radius of Rp ~ 0.35A, 

Ff = Ek=i AFp{Rl/Rp)eiiRp - i?r,)/A)Rr,. Here 9 is 
the Heavisidc step function, IV^ is the location of pinning 



site fc, 



andRf, = (R, 



The ex- 



A{F^k + Fy-y) represents the Lorentz 



Ri 

ternal force Ff = 

force from an applied current. The vortex velocities Vx 
(Nv)-^ ^ • dRJdt and Vy = (N^)-^ ^^I", y • dRJdt 
in the two directions represent the resulting voltages. In 
this work, all driving forces are small enough that the vor- 
tices at the pinning sites remain pinned. The initial vortex 
positions are found by simulated annealing in a procedure 
similar to that used previously [9]. After the vortex po- 
sitions are initialized, we apply a fixed force F^ in the 
positive y-dircction while slowly increasing a force F!^ in 
the positive x direction. 

We first consider the vortices under matching conditions 
B = B^, where the number of vortices equals the number 
of pinning sites in the original square pinning array. In 
the absence of a drive, the vortices form a square lattice 
with the interstitial vortices sitting at the locations of the 
pinning sites that were removed. We illustrate the vortex 
trajectories for several values of FJjr' at fixed F^ = 0.2 
in fig. 1, and we plot versus Fjj^ in fig. 2(a), dV^/dF^ 
versus F^ in fig. 2(b), and Vy versus F^ in fig. 2(c) for the 
same system. The points marked a-e in fig. 2 correspond 
to the values of F^ illustrated in fig. l(a-e). 

At F^ = 0, the vortices in the y channel are flowing 
while the vortices in the x channel are pinned. The vor- 
tices in the x channel do not depin until F^ reaches a crit- 
ical value F^ = 0.037. For Ff < 0.037 the vortices in the 
X channel are pinned but those in the y channel are mov- 
ing, as indicated by the fact that T4 = over this regime 
in fig. 2(a) while Vy has a finite value in fig. 2(c). Fig- 
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Fig. 2: (a) 14 vs , (b) dV^./dF^r' vs F^ , and (c) Vy vs F^ 
for the system in fig. 1 with F^ — 0.2. The points marked 
a-e correspond to the flow regimes illustrated in fig. l(a-e). 
Upper left inset in (a); depinning threshold F^ for motion in 
the s-direction for varied values of F^ for the same system. 
Lower right inset in (a): a blowup of 14 vs F^ from the main 
panel highlighting the decrease in Vx with increasing F^ as 
the system exits the F^/F^ — 1.0 locked phase. In (b) sev- 
eral regions of negative differential conductivity appear with 
dVx/dF^ < 0. 



ure 1(a) illustrates the vortex trajectories for = 0.025 
where the flow is confined to the y direction. At = 0.4, 
fig. 1(b) shows that vortices in both channels arc moving 
and that the trajectories arc disordered at the channel in- 
tersection. 

In fig. 1(c) at F^ = 0.106, the vortex flow at the inter- 
section is ordered or coherent and consists of a mixing flow 
where vortices approaching the intersection from the left 
turn upwards into the y channel, while a portion of the vor- 
tices approaching the intersection from below turn right 
into the x channel. The net flow at the intersection con- 
sists of two vortices exiting in the positive y direction for 
every one vortex exiting in the positive x direction. The 
ordered flow regime correlates with a steplike transport 
anomaly centered around F^ = 0.1 at the point marked c 
in fig. 2(a). Two peaks in dVx/dF^ in fig. 2(b) surround 
this region, while a smooth rounded peak appears in Vy in 
fig. 2(c). This shows that the vortices remain locked in the 
ordered flow phase illustrated in fig. 1(c) for a range of F^P. 
In general, ordered flow phases are associated with all of 
the transport anomalies in flg. 2. The phases are centered 
at rational ratios of the external drives F^ /F^ = n/m 
with n,TO integer, such as at F^ /F^ — 1/4, 1/2, or 1.0. 
Figure 1(c) shows the F^ / F^ = 1/2 state. 

In fig. 2, anomalies at other rational ratios of Fy /F^ 
are visible, with the most pronounced steplike features 
occurring for cases where n < A and to < 4. This type 
of behavior is very similar to the mode locking found in 
systems which contain an intrinsic oscillating frequency 
that locks to an external ac drive [13]. Another intrigu- 
ing feature is the fact that near the drive at which the 
system exits a locked phase, such as near the ends of the 
Fy /F^ = 1/2 and 1.0 steps, Vx decreases with increasing 
F^ . This is highlighted in the lower right inset of fig. 2(a) 
for the end of the F° /F^ = 1.0 locked phase. The ve- 
locity decrease produces negative differential conductivity 
where dV^/dF^ < 0, as shown in fig. 2(b). This indicates 
that in order for the system to remain locked as it ap- 
proaches the upper end of the locked phase, the vortices 
must move more slowly in the x direction as FJl? increases. 

The anomalies we find in the velocity force curves are 
similar to those observed for vortices or colloids driven 
over periodic pinning arrays when the particle velocity 
vector locks to symmetry directions of the pinning ar- 
ray [16-18]. In these systems, the locking produces clear 
steps and occurs even in the limit of an isolated driven 
particle [18]. The phase locking considered here for the 
crossed channel geometry results from a different mecha- 
nism and docs not occur in the single particle limit, as we 
describe below in fig. 3. In addition, there was no negative 
differential conductivity associated with the phase locking 
in rcfs. [16-18]. 

Away from the steplike velocity anomalies, the vor- 
tex trajectories are disordered, as shown in fig. 1(d) for 
= 0.1105. Figure 1(e) illustrates the ordered flow 
phase at F,^ = 0.2, which corresponds to the pronounced 
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Fig. 3: (a) Vx vs for the same system in fig. 2(a) but with 
Fjf = 0.1. The initial pinned phase and the F°/F° = 1.0 
state are enhanced, (b) 14 (dark fine) and Vy (light line) vs 
F^ at _Fjf — 0.7 with only a single interstitial vortex in the 
channels. All the anomalies seen in fig. 2 are lost and there is 
only a sharp transition from y direction motion to x direction 
motion as a function of FF . 



/F^ = 1.0 anomaly in fig. 2. Also associated with this 
ordered flow is the dip in Vx highlighted in the lower inset 
of fig. 2(a), the arctangent-like feature and large negative 
differential conductivity in dVx/dF^ in fig. 2(b), and the 
smooth feature in Vy at point e in fig. 2(c). Figure 1(f) 
shows the disordered vortex trajectories at F.^ ~ 0.224, 
just outside of the Fy /F^ = 1.0 state. 

The general features we find in the transport character- 
istics are independent of the rate at which F^ is increased, 
and we find similar features if we vary F^ at fixed F^ . 
In the upper left inset of fig. 2(a) we show the depinning 
force F!: for motion in the x channel vs F.f . As F.f in- 
creases from zero, F!^ initially decreases and then begins 
to increase. The increase occurs when the y direction driv- 
ing force shifts the vortices in the x channel closer to the 
occupied pinning sites along the upper side of the chan- 
nel, increasing the coupling between the interstitial and 
pinned vortices and thus increasing the depinning force for 
the interstitial vortices. The widths of the velocity anoma- 
lies are enhanced when the depinning threshold increases. 
This is illustrated in fig. 3(a) where we plot Vx versus F^ 
for the same system in fig. 2(a) but with F^ = 0.1. Both 
the pinned phase and the F^ /F^ = 1.0 state are broader 
at Fjf = 0.1 than at F^ = 0.2. 

We next examine the case of a single interstitial vortex 
driven through the channels. We remove all but one of 
the interstitial vortices and initialize this single intersti- 
tial vortex in the y channel. In fig. 3(b) we plot Vx and Vy 
versus F^ at F^ = 0.7 showing that in the single particle 
limit, all the transport anomalies seen in fig. 2 are lost and 
there is simply a transition from pure y direction motion 
to pure X direction motion when F^ increases. This is 
indicated by the sharp transition to Vy = 0.0 which coin- 
cides with the onset of a finite value of Vx ■ We term this 
sharp transition a vortex gating effect. 

The loss of the phase locking for a single interstitial vor- 
tex indicates that collective interactions are responsible for 
the phase locking in the crossed channel geometry. The 
discreetness of the vortices arriving at the intersection of 
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Fig. 4: Vortex positions (filled dots), pinning site locations 
(open circles), and trajectories (lines) for the same system in 
fig. 1 but with a higher vortex density of B/B'^ — 1.067. Here 
the vortices in the channel undergo a buckling transition to a 
zig-zag pattern, (a) Ordered fiow at the F^ /F^ = 1.0 state, 
(b) Disordered fiow at F° jF^ = 0.92. 



the channels permits only an integer number of vortices to 
move in each of the two directions. The ordered trajecto- 
ries such as those seen in fig. l(c,e) distort somewhat but 
stay ordered as F^ increases, allowing the interstitial vor- 
tices to remain in the same phase-locked state over a wider 
range of Fj^ . The distortion of the vortex trajectories is 
also responsible for the negative differential conductivity 
and the changes in Vy associated with the coherent flow 
phases. In the case of vortices or colloids moving over peri- 
odic pinning sites in Refs. [16, 18], the ordered flow phases 
consisted of one-dimensional paths oriented along different 
pinning lattice symmetry angles. These paths cannot dis- 
tort and so there was no differential negative conductivity 
in these systems. 

The phase locking illustrated in figs. 1 and 2 occurs 
for a range of fields and is most prominent for 0.925 < 
BjB^ < 1.1. Within this field range, the interstitial vor- 
tices are confined in the channels. In fig. 4(a,b) we plot 
the vortex trajectories for a sample with B/B^ = 1.067. 
In this case, the vortices in the channel undergo a buckling 
transition and form a zig-zag pattern, fig. 4(a) illustrates 
the F^j^ /F^ = 1.0 state where an ordered flow pattern 
occurs similar to that seen in fig. 1(e). Here there are two 
staggered moving rows of vortices in each channel. For 
Py /Px = 0-92, shown in fig. 4(b), the vortex trajecto- 
ries are disordered and there is no steplike anomaly in the 
transport curves at this point. For B/B^ > 1.1, some 
interstitial vortices begin to penetrate the bulk pinned re- 
gion away from the channels. These interstitial vortices 
are generally more strongly pinned than the vortices in 
the channels; however, for intermediate ranges of F^ and 
F^ they can become mobile and add extra features to the 
velocity-force curves. These effects are beyond the scope 
of this work. 

In summary, we studied a square periodic pinning array 
with crossed easy flow channels for vortex motion formed 
by the removal of one row and one column of pinning. We 
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show that when a fixed drive is applied in one direction and 
a slowly increasing drive is applied in the perpendicular 
direction, a series of phases with coherent vortex flow dy- 
namics occurs which are centered at rational ratios of the 
two drives. These coherent phases appear as anomalies in 
the transport curves and are also associated with negative 
differential conductivity. The phase locking anomalies are 
distinct from those observed for the velocity vector locking 
of particles driven at different angles over periodic pinning 
arrays, which appear even in the single particle case. For 
the crossed channel geometry, the phase locking occurs 
due to the discrete nature of the particles. In the single 
particle limit, the transport anomalies vanish and are re- 
placed by a vortex gating phenomena in which the motion 
switches abruptly from one channel to the other. We ex- 
pect that similar effects can occur for particle transport 
in other systems with quasi-one-dimcnsional crossed ge- 
ometries, such as colloidal systems, microfluids, or crossed 
one-dimensional Wigner crystal systems. 

^ ^ ^ 

This work was carried out under the auspices of the 
NNSA of the U.S. DoE at LANL under Contract No. DE- 
AC52-06NA25396. 
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